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Testing universality and the fractional power prescription for the 
staggered fermion determinant 

David H. Adams a * 

a Instituut-Lorentz for Theoretical Physics, Leiden University, 
Niels Bohrweg 2, NL-2333 CA Leiden, The Netherlands 

In Ref. U expressions for the continuous Euclidean time limits of various lattice fermion determinants were 
derived and compared in order to test universality expectations in Lattice QCD. Here we review that work 
with emphasis on its relevance for assessing the fractional power prescription for the determinant in dynamical 
staggered fermion simulations. Some new supplementary material is presented; in particular the status of the 
"universality anomaly" in the determinant relations is clarified: it is shown to be gauge field-independent and 
therefore physically inconsequential. 



1. Introduction and background 

During the last decade there have been major 
developments in lattice gauge theory. At the con- 
ceptual level, the long-standing problem of how 
to formulate chiral symmetry on the lattice has 
been solved. The explicit solution is provided by 
the Overlap fermion formulation which led 
to the Overlap Dirac operator for LQCD [3] - 
this has an exact chiral symmetry encapsulated 
in the Ginsparg- Wilson relation 0] 75D + D^ 5 = 
aDj 5 D, and is free |S] of the problems caused by 
lack of chiral symmetry in previous formulations. 
Implementing Overlap fermions in numerical sim- 
ulations is a very difficult challenge though, and 
it seems that substantial advances in algorithms 
and computer power will be needed before realis- 
tic simulations with dynamical Overlap fermions 
are possible. 

In the meantime work continues on traditional 
fermion formulations. A major practical develop- 
ment has been algorithm advances for improved 
staggered fermions which have made possible dy- 
namical fermion simulations on realistically large 
lattices. An impressive increase in agreement 
with experimental values for various hadronic pa- 
rameters has been obtained compared to previ- 
ous quenched simulations [S] (see, e.g., |7I8| for 
reviews and [5] for latest results). However, as 
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discussed previously by various authors [1018111) . 
there are problematic conceptual/theoretical is- 
sues for dynamical staggered fermions which raise 
the question of whether this is a hrst-principles 
approach to QCD or whether it should instead 
be regarded as some kind of effective model. 

The staggered formulation is a lattice theory 
for 4 degenerate continuum fermion flavours, and 
to represent the fermion determinant of a sin- 
gle dynamical (sea) quark flavour one takes the 
fourth root of the staggered fermion determinant. 
This prescription is based on the expectation that 
in the continuum limit the staggered fermion de- 
terminant should factor into the product of 4 
copies of the fermion determinant for a single 
quark flavour. There are two theoretical issues 
here: firstly, whether the expected factorisation 
actually occurs, and secondly, assuming it does 
occur, whether the prescription can be fitted into 
the framework of local quantum field theory at 
finite lattice spacing, i.e. whether there exists 
a (exponentially-)local lattice Dirac operator D 
such that det(D staggered) 1 ^ = detD. Here we 
will focus on the factorisation issue. For discus- 
sion of the locality issue see [1018111) along with 
§IX-D.7 of [§| and the references mentioned there. 

It is well-known that terms in lattice actions 
which formally vanish for a — > can have resid- 
ual effects which remain in the continuum limit. 
Familiar examples include non-vanishing contri- 
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butions to Feynman diagrams from "irrelevant" 
interaction terms in lattice perturbation theory 
(see, e.g., ^21) an d the essential role of the Wilson 
term in reproducing the axial anomaly in LQCD 
with Wilson fermions I n light of this, factori- 
sation of the staggered fermion determinant in the 
continuum limit cannot be taken for granted - it 
could conceivably be spoiled by residual effects 
from flavour-changing interactions. It is there- 
fore important to investigate, to the extent that 
it is possible and in as direct a way as possible, 
whether the factorisation actually occurs. 

Recall that the expansion of the logarithm 
of the fermion determinant in powers of the 
gauge field can be expressed in terms of the 
one (fermion) loop gluonic n-point functions. 
Therefore, the factorisation issue for the stag- 
gered fermion determinant is intimately tied to 
the question of whether perturbative LQCD with 
staggered fermions can be renormalised in a way 
that is consistent with 4 flavour QCD. The stud- 
ies of perturbative LQCD with staggered fermions 
carried out to date indicate that this is the case, 2 
and this can be taken as indirect evidence for fac- 
torisation of the staggered fermion determinant 
in the continuum limit (although it says nothing 
about possible non-perturbative problems for the 
factorisation) . It should be noted however that so 
far renormalisation of perturbative LQCD with 
staggered fermions to all orders in the gauge cou- 
pling has not been demonstrated. (Some argu- 
ments were sketched in but these are not a 
rigorous proof.) This is in contrast to the situ- 
ation for Wilson fermions where renormalisabil- 
ity was demonstrated some time ago by Reisz 
[TC] based on his power-counting theorem [T7| . 
As mentioned in J0|, the power-counting theo- 
rem, as it stands, does not apply to staggered 
fermions. This reflects the additional complica- 
tions due to flavour-changing interactions in stag- 
gered fermion perturbation theory. In light of all 
this, investigation of the factorisation issue for the 
staggered fermion determinant is also of interest 
as an indirect check on the possibility of renormal- 

2 Early indications of this came from evaluations of the 
fermion loop correction to the gluonic propagator 1141 . and 
the fermion self-energy 1151 : see 1121 for a summary of 
subsequent work. 



ising LQCD with staggered fermions to all orders 
in a way that is consistent with 4 flavour QCD. 

The factorisation issue for the staggered 
fermion determinant is closely tied to a more gen- 
eral universality issue: Is LQCD with a staggered 
fermion in the same universality class as, say, 
LQCD with 4 flavours of Wilson fermions? This is 
certainly expected to be the case but should nev- 
ertheless be tested wherever possible. One way to 
try to get insight into this, while at the same time 
investigating the factorisation issue, would be to 
compare the staggered fermion determinant with 
the fourth power of the Wilson fermion determi- 
nant in the continuum limit. Ideally this should 
be done analytically (so as to be able to iden- 
tify and discard physically inconsequential factors 
which may be present in the determinants and 
which may either diverge or vanish for a — > 0); 
however this appears impossible with currently 
known techniques. A simplified version of this 
test of universality is possible though pQ and we 
review it in the following. When supplemented 
with some additional results presented here, it 
allows us to analytically study the factorisation 
issue for a "partially staggered" fermion determi- 
nant in a setting where only a partial continuum 
limit is required. From the results of U (and the 
clarification of the "universality anomaly" which 
we give here) we find that factorisation of the par- 
tially staggered fermion determinant does indeed 
occur in this limit. 

2. A simplified framework for testing uni- 
versality 

On a finite volume Euclidean spacetime lat- 
tice with temporal lattice spacing a and spacial 
lattice spacing a' we consider the fermion ac- 
tion S fermion = a(fl') 3 Z)(x,r) ^( X > T)D^1p(x, t) 

where 

D (r) = i 74 V 4 + ^A 4 + D space + m (I) 
n — _L j_ r' a 

space — a'' space 1 2a' space 
with. space — ^2(7—1.2 3 Tc^cr j ^space 

E CT =i, 2 ,3 A ° > V n ip(x) = ^{U^x^ix+bfiyU^x- 
bp)^(x-bp)] , A u ^(x) = 2tp(x)-U n (x)ip(x+bp)- 
U a (x — bp)^ 1 ip(x — bp,) , x— (x,t) and b = a or 
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a' as appropriate. For r = r' ^ this is the 
Wilson fermion action, while for r = r' = it is 
the naive action. In the case r — 0, r' ^ it is 
a partially naive — partially Wilson action, with 
fermion doubling on the Euclidean time axis due 
to the absence of the temporal part of the Wil- 
son term. Universality can therefore be tested 
by comparing this lattice theory with the theory 
for two flavours of Wilson (r ^ 0) fcrmions: ac- 
cording to universality these should coincide in 
the continuous Euclidean time limit a — * 0. We 
do not need to take the spacial continuum limit 
since the spacial parts of the actions of the two 
theories are the same. Thus we have a simplified 
framework for testing universality which only re- 
quires taking the continuum limit of one of the 
spacetime coordinates (chosen here to be time). 
The key advantage is that, as we will see, this is 
accessible to direct analytic investigation. 

To connect this universality test with the fac- 
torisation issue we point out that the r = theory 
has a partially staggered interpretation as follows. 
Introduce the staggered 'flavour' fields 

V»i(x,r) = ^(x, 2t) , ip2(x, t) = ^(x, 2t + a) 

living on the partially blocked lattice with tempo- 
ral spacing 2a. In the free field case £>( r=0 ) then 
takes the form 



n(0) ( V'i 



i74 



(»+ 



(£) 



where dj (d^) denotes the forward (backward) 
time difference operator. Now, making a basis 
transformation in spinor®flavour space specified 

by = (_ 7 ^ 4 \) we obtain 



+D S 



'^1 J 2(2a) 

y >/ ( ;i,T (2) 

where = h{9± + dj). This clearly has the 
structure of a 'partially staggered' lattice Dirac 
operator: the 'time part' has the same form 
as the usual free field staggered Dirac opera- 
tor in the flavour field representation |18j . In 
particular, A4 comes with a non-trivial flavour 
matrix, and in the gauged theory this gives rise 
to flavour-changing interactions. Thus we have a 



simplified setting for investigating the question of 
whether residual effects from flavour-changing in- 
teractions spoil the factorisation of the staggered 
fermion determinant in the continuum limit. In 
the present setting, factorisation of the partially 
staggered fermion determinant in the continuous 
time limit corresponds to the following 



Universality expectation: 

lim aletD( r =V = ( i im detD^V 

a->0 V a^O 



(3) 



(up to physically inconsequential factors 
(p.i.f.'s)). 

This is something which can be checked analyti- 
cally P as we describe in the following sections. 

3. Results 

Set j3 :=length of the time axis (which we 
hold fixed when taking the continuous time limit 
a — > 0); f3 — aNp where Np is the number of 
sites along the time axis; set N :=dimension of 
the quantum-mechanical Hilbert space of spinor 
fields {^(x)} living only on the spacial lattice, 
and take U^x, r) to be the lattice transcript of 
the 4-component (x, r) of a smooth continuum 
gauge field. In this setup the following was shown 
in [T] (full details in H2j): 
Result of direct calculation: 



detD {l 



detD<® 



i f 3 TrMMdr 

: 2 Jo ' 



, 2a 



,NN B 



det(l 



det(l-V(fl)) 
(4) 

V(/3)) 2 (5) 



with the ingredients defined below. The gauge 
fields are assumed to satisfy periodic time bound- 
ary condition, and for simplicity we have stated 
the results (4)-(5) for the case where the fermion 
fields are also time-periodic. (The treatment in 
P covers the case of general time b.c. ^( x : r) = 
e -Q/3 V( x jO). p or a — fx + in/p this corresponds 

to QCD at finite temperature and chemical 
potential /1.) In (4), 

M(t) := ^A space {T)+m (6) 

is a linear map on the QM Hilbert space 
{^(x)}, with A spoce (r)V'(x) given simply by re- 
placing ^>(x, r) with V( x ) in the definition of 
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A spaced '( x i r )- The V(/3) in (4)-(5) is also an 
operator on the QM Hilbert space; it is the spe- 
cialisation to t = (3 of the operator V(r) defined 
Vr G R as follows. We regard V^X; 7 ") as a func- 
tion \&(t) taking values in the QM Hilbert space, 
and introduce the continuous time — lattice space 
Dirac operator 

D = l4 (-§ f +A 4 (T))+D space (T)+m (7) 

After extending the gauge field from r G [0, (3} to 
all r € R, periodic under r — > r+/3, this operator 
acts on QM Hilbert space- valued functions ^(r) 
defined for all r G R. Solutions to D*(r) = 
(t G R, no periodicity requirement on \&(t)) are 
then specified by an initial value \P(0), and V(r) is 
defined to be the evolution operator which deter- 
mines the full solution from the initial value, i.e. 
\fr(r) = V(r)*(0). This completes the descrip- 
tion of the ingredients in (4)— (5). We outline the 
derivations of (4)-(5) in the next section. 

The results (4)-(5) allow us to check to what 
extent the universality expectation (3) is satis- 
fied. The appearance of det(l — V(/3)) in (4) ver- 
sus det(l — V(/3)) 2 in (5) is clearly in accordance 
with this expectation. The divergent factors in 
(4)-(5) are also in accordance with it: given that 
the divergent factor in (4) is (\/a) NNf3 , and in 
light of the partially staggered interpretation of 
the r = theory, the divergent factor in (5) should 
be ((l/a)^) 2 where a = 2a and N = Np/2. 
This is precisely the divergent factor (l/2a) NN/3 
appearing in (5). However, (4)-(5) reveal a po- 
tential breakdown of the universality expectation 

due to the factor e-fo TrM ( T ) dr appearing in (4) 
which has no counterpart in (5). Thus it is im- 
portant to clarify the physical significance, or lack 
thereof, of this "universality anomaly" . Recalling 
(6) we see that, up to a p.i.f., the anomaly factor 
is 

exp{± / Tr(£rA space ( T ))dT} (8) 
Jo 

This involves the spacial link variables and can- 
not therefore be immediately dismissed as a p.i.f. 
One might argue that it becomes a p.i.f. when 
the spacial continuum limit a' —> is taken since 
-^■A space formally vanishes in this limit. This is 



a delicate issue though, since Tr(-^A space ) ac- 
tually diverges in this limit (the largest eigen- 
value of -^rA space is ~ ■%)• In H] the status of 
this anomaly factor was left unresolved. We now 
clarify the situation, showing that the anomaly 
is in fact physically inconsequential. This is 
done by showing that TrA space (r) , and hence the 
anomaly itself, are independent of the gauge field, 
i.e. do not depend on the spacial link variables 
entering in A spoce (r). To see this, recall that 
A space (Y) = E CT =i,2,3 Aff(r) with A ^( T Mx) = 
2ip(x) - U a (x, T)V>(x + a'<7) - C/ .(x-o'o-) -1 ^(x- 
a'a). Evaluating TrA a (r) in a spacial plane wave 
basis "0p( x ) ~ eipx we S e t (ignoring overall nor- 
malisation factors) 

Tr A CT (r)=^(^ p ,A^p> = 
p 

P x 

After interchanging the sums over p and x the 
terms involving the link variables are seen to van- 
ish since Y. Pa eV " a ' =S(a')=0. 

Thus we have established that the universality 
expectation (3) does indeed hold up to physically 
inconsequential factors, at least when r = 1 in the 
right-hand side of (3). This is at the same time 
a demonstration of the factorisability of the par- 
tially staggered fermion determinant in the par- 
tial continuum limit a — > 0. For technical reasons 
mentioned below we have only been able to eval- 
uate the a — * limit of detD^ r ' in the cases r = 0, 
Eq.(5), and r=l, Eq.(4), but not yet for general 
r ^ 0. 

4. Derivation of the results an outline 

The time-periodic lattice spinor fields r) 
on which acts are identified with QM Hilbert 
space- valued functions ^(r) living on the tempo- 
ral lattice sites, and these can in turn be rep- 
resented by vectors * = (*(0), . . . , $(Np - 1)) 
where ~&(k) := ^(ak). Then is represented 
by 

Z)W^(fc) = 

d£\ (fc)*(fc- 1) + 4 r) (*o*(fc) + 4 r) (fc)*(fc+i) 
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space(k) + m with Ui{k) := U^ak) and 
,,,ace(k) := D space (ak). (Here U 4 {ak)ip(yL) := 



where the dj (k) 's are linear operators on the QM 
Hilbert space given by d^\k) = ^(74 — rW^k) 
d -l( k ) = -^(l4 + r)U 4 (k-li 

b 

U4(x,ak)ip(x) etc.) After writing Z)( r ' as an 
N/3 x N/3 matrix, its determinant can be straight- 
forwardly evaluated via the method of [201 ■ The 
cases r = ±l and r =^ ±1 require separate treat- 
ments due to the fact that d±\(k) is invertible 
when r / il but not when r = ±1. For sim- 
plicity we restrict to r > (the r < case is 
analogous), then the results of the determinant 
calculations are (with Np even in the r ^ 1 case) 

nnni: 

{ i^pif^det((l (9) 
detC« = (i)^ x(M) <fet(l - V{Np)) (10) 

where x(M) := Jlfc^o" 1 rfe *(! + aAf (0A;)) 1 / 2 and 
W( r '(n) and V(fc) are the lattice evolution opera- 
tors in the r ^ 1 and r = l cases, respectively: 

D^)*(*)=o <* (4 ( :? 1) )=^ ) w(is) 

£)W*(fe) = ^ *(fc) = V(Jfe)#(0) 

(A: € Z, no periodicity requirement on $(k)). 
Note that in the r ^ 1 case two initial values are 
needed to determine the solution, hence W r \n) 
is a 2 x 2 matrix whose entries are linear operators 
on the QM Hilbert space, whereas in the r = 1 case 
it turns out that only one initial value is needed. 

The next step is to show that the lattice time 
evolution operators converge to appropriate con- 
tinuous time evolution operators in the a — > 
limit. The idea is to write D^ r ' in the form 



£>W = L 1 (fc)ia + L (fc) 



(11) 



or some suitable variant thereof. Here d is the for- 
ward or backward finite difference operator and 
Lj(k) (j — 1,2) are linear operators on the QM 
Hilbert space, parameterised by k £ Z, periodic 
under k — > k+N/3 , such that 



for some operators Lj(r) parameterised by con- 
tinuous time variable tgR with periodicity un- 
der r -> t+/3. Then the solutions to Z)M#(fc) = 
approximate the solutions to D^>(t) = where 



= £i(t)£+L O (t) 



(13) 



Consequently the lattice evolution operator U(k) 
for solutions to = approximates the continu- 
ous time evolution operator U(t) for solutions to 

D* = 0. (Explicitly, U(t) = Te~ So M^M*)* 
where T — i-ordering.) In particular one has 
the convergence theorem: lim _>o tl(Np) = U(f3). 
This type of convergence result is well-known in 
the mathematics literature [21] • It is the key to 
deriving the results (4)-(5) from the determinant 
formulae (9)-(10). 

We have so far only been able to write 
in the form (11) in the r = and r = 1 cases. 
The explicit expressions in these cases are given 
below. The problem of evaluating lim _»o detD^ 
for general values of r remains for future work; 
new techniques beyond those described here may 
be required for this. 

r = 1 case : Decomposing \&= ^ , 74^± = 
±i&± , we have 



Li(k) 
L (k) 



Li \ ( a a°+) + 
74 V i(Ui(k)-i] 



-D space (k) + m 



The Lj(k)'s are periodic under k — > k + N/3 
and satisfy (12) with L\{t) — 74 and Lq(t) = 
74^4(7-) + D(t) + m, hence the continuous time- 
lattice space operator (13) is precisely the opera- 
tor D introduced earlier in (7). The convergence 
theorem then gives lim a _»o V(Np) — V(/3). The 
claimed result (4) now follows from (10) after not- 
ing that lim a ^ X (M) = e* So Tr M(r) dT [T]. 



Lj(k) = L J {ak) + 0{a) 



(12) 



/ *i(n)\ = / *(2n) \ 
V*2(n)/ \^{2n+l)) 

- f 1 ( d+ 



f=)(r=0) _ 



we have 



in 
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L Q (n) = 

/ D space (2n)+m 74 £ (£/ 4 (2n+l) -U A (2m)" 1 ) 

(f>4(2ri)-f> 4 (2n-l)- 1 ) £> space (2n)+m 

The Lj(n)'s in this case are periodic under n — > 
n+Nfj/2. Furthermore, £j(n) = L 3 (2an) + 0(a) 
with the corresponding continuous time — lattice 
space operator (13) given by 

~ = / D space { T ) + m 74 (^ + A 4 (r)) \ 
V 74(^7 + ^4(r)) D space {T) + m J 

This can be rewritten as£) = -1 ^ ^0 with 

D as in (7) and as defined above Eq.(2). The 
convergence theorem can now be applied to get 
lim 5 ^o U {r=0) (Np) = 0(T V ( P) )&-\ where 

a = 2a and Np = N/3/2. This together with (9) 
gives the claimed result (5). 

Detailed rigorous proofs of the convergence the- 
orem in the specific cases where it has been used 
above are given in [T5|. 

5. Free field case 

In the free field case the evolution operator 
V(j3) has a simple expression: V(j3) = er^ H 
where H is the QM Hamiltonian, i.e. H = 

l^Dgpace — liOtf space + M) , M — ^p-A space + 

m. The energy eigenvalues of H are ±E(p) = 
y/p 2 + M(p) (2- fold degenerate), and it follows 
that 

det(l-Vm = n[(l-0(l- e -^)] 2 
p 

= l[[e^(l-e-n 2 } 2 
p 

Substituting this into (4)-(5) we obtain expres- 
sions for the a — > limits of detD^ and 
detD {Q) in terms of the £(p)'s. These expres- 
sions can be re-derived starting from debDv> — 
{product of eigenvalues] , using Matsubara fre- 
quency summation techniques |19| . This provides 
a crosscheck on the correctness of the results. 



6. Ambiguous continuous time limit of the 
partially staggered fermion determinant 

In obtaining the formula (9) for detD^ r=0 ^ the 
number Np of lattice sites along the time axis 
was assumed to be even, and it was under this 
restriction that the a — > limit result (5) was 
obtained. One might expect that the a — > limit 
of detD^ with odd Np is the same. However, 
this turns out not to be the case. For odd Np one 
finds ^21, instead of (9), 

detD^ = 

(^) NN ?det[(l \)-{\l)U^{Np/2)] 

where the precise meaning of U i ~ r) {Np/2) in the 
odd Np case is given in ^] and we simply men- 
tion here that U^ 0) (Np/2) has the same a — > 
limit as in the even Np case. It then follows by a 
straightforward calculation ^U] that in this case 

defcDCO) det(l-V{p))det{\+V{p)) 

Comparing with (5) we see that one of the fac- 
tors det(l — V(/3)) in (5) is replaced here by 
det(l + V(/3)). Thus the continuous time limit 
of the partially staggered fermion determinant 
detD^ is ambi guous: different answers are ob- 
tained depending on whether Np is restricted to 
be even or odd. 

This ambiguity has a natural explanation 
though: it reflects the fact that even Np and odd 
Np in the r = lattice theory correspond to dif- 
ferent time boundary conditions in the underly- 
ing continuous time theory. This is seen from the 
partially staggered 'flavour field' interpretation as 
follows. For even Np the flavour fields are time- 
periodic just like the original field, but for odd 
Np we have 

*i(0) = #(0) = $(Np) = M{Np - l)/2) 
* a (0) = (1) = 4f(Np + 1) = &i{(Np + l)/2) 

(Ho!) - (^SKg). which in the 
a^01imitbecomes(^)) = (° 
Thus the b.c.'s in the underlying continuous time 
theory are "twisted" in this case. 
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7. Conclusions 

We have reviewed a test of universality where 
a "partially naive" lattice fermion theory with 
fermion doubling on the Euclidean time axis 
is compared to the theory for two degenerate 
flavours of Wilson fermions £Q. The continu- 
ous Euclidean time limits of the fermion determi- 
nants were evaluated and compared (with Wil- 
son parameter r = 1 for the Wilson fermion). 
From the expressions obtained, the partially naive 
fermion determinant was seen to coincide with the 
2 flavour Wilson fermion determinant (i.e. the 
square of the single flavour Wilson fermion deter- 
minant) up to physically inconsequential factors, 
thus confirming the universality expectation. In 
connection with this the "universality anomaly" 
found in was shown to be physically inconse- 
quential. 

Furthermore, the partially naive theory was 
shown here to have a "partially staggered" in- 
terpretation, thus our result can be interpreted 
as demonstrating factorisation of a partially stag- 
gered fermion determinant on the partial contin- 
uum limit a — > 0. This is a first piece of ana- 
lytic evidence suggesting that factorisation of the 
complete staggered fermion determinant occurs in 
the full continuum limit (a problem with factori- 
sation in the complete case could very well also 
show up as a problem in the partially staggered 
case). In light of the discussion in §1, our result 
can also be interpreted as providing a hint that 
renormalisation of perturbative LQCD with (at 
least partially) staggered fermions to all orders, 
consistent with 4 flavour QCD (or 2 flavour QCD 
in our partially staggered case), is possible. 3 

I would like to thank the organizers of Lattice 
2004 for a very enjoyable and stimulating confer- 
ence, and for the opportunity to present this work 
in the plenary session. I also thank the NCTS 
at Taiwan National University for hospitality and 
support during visits where some of this work was 
done, and Pierre van Baal for feedback on the 
manuscript. 



3 We remark that renormalisation of a variant of partially 
staggered fermions in the context of a chiral sigma model 
(an effective theory; no gauge fields present) has been car- 
ried out in 1221 . 
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